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Abstract 

We discuss qualitative behavior of the SU(N) gauge beta functions in QCD with many massless 
flavors. Non-perturbative beta functions can be obtained by extracting the renormalized tra- 
jectories in the exact renormalization group framework. We examine the renormalization group 
equations for the general four-fermi couplings as well as the gauge coupling obtained in a simple 
approximation scheme. It is shown that the gauge beta function possesses not only an IR but also 
a UV fixed point in the conformal window. These fixed points merge with each other and disappear 
at the edge of the conformal window. The scaling dimensions of the quark mass operator at the 
fixed points are also shown. 
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I. INTRODUCTION 



Realization of scale invariant quantum field theories in four dimensions is a very attracting 
subject not only as a theoretical interest but also for the gauge hierarchy problem of the 
standard model. It has been known for sometime that non-abelian gauge theories with 
appropriate numbers of massless fermions can be scale invariant. Explicitly, the IR attractive 
fixed point, called the Banks-Zaks (BZ) fixed point has been shown to realize in the beta 
function of an SU(N C ) gauge theory, as long as the flavor number Nf is less than but close to 
lliVc/2. Then the fixed point appears in the weak coupling region and perturbative analysis 
is reliable. As the flavor number is reduced, the fixed point coupling becomes strong and 
enters into the non-perturbative region eventually. 

The non-perturbative dynamics of QCD with many flavors has been studies by analyzing 
the so-called ladder Dyson-Schwinger (DS) equations It has been known that the 

chiral symmetry is spontaneously broken for the gauge couplings stronger than a certain 
critical value. This implies that the IR fixed point does not exist beyond the critical gauge 
coupling, since scale invariance is lost by the dynamical scale generation. The chiral sym- 



metry breaking has been also examined [5] by means of the Wilson renormalization group 
jfj], or more explicitly by solving the exact (functional) renormalization group (ERG) equa- 
tions Jfj-lO]. The analyses using the ERG equations have been also performed for the chiral 



dynamics in the conformal window 
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Thus it is expected that the lower boundary of the conformal window N^ r is determined 
by spontaneous breakdown of the chiral symmetry. The DS equations in the large N and 
adder approximation leads to Nj r = 4iV c , when combined with the two-loop beta function 
2)]. A recent ERG analysis combined with the four-loop gauge beta function also leads to 
N^ r = lO.Ojj^y for the SU(3) case ll|. Lattice Monte-Cairo simulations have been also 



performed 14J and the results indicate existence of the conformal phase and the critical 
flavor number. The recent simulations have provided evidence that 8 < Nf" < 12 for the 
SU(3) QCD, even though the case Nf = 12 is controversial [la, Il6 |. 

In this paper we consider how the fixed point disappears from the gauge beta function 
when the flavor number goes out of the conformal window. The IR fixed point in the two- 
loop perturbative beta function remains even below Nf and seems to move away to infinity 
at a certain Nf lower than Nf"- Meanwhile chiral symmetry breaking at the strong coupling 
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region prohibits existence of the fixed point. Then it seems that the beta functions in the 
conformal window cannot be continued smoothly to the beta function with a fewer flavor 
number 34]. 



The unique possibility is that a UV fixed point exists as well in the conformal window 
and merges with the IR fixed point at the edge of the conformal window. Then the IR fixed 
point can disappear before reaching the critical gauge coupling for chiral symmetry breaking. 



Such mechanism has been discussed in Ref. [17] and showed that the BKT scaling, or the 
so-called Miransky scaling fl Q for the conformal gauge theories, realizes. In this paper 
we will show, although in a certain scheme of approximation, that the gauge beta function 
possesses a UV fixed point and fixed point merger occurs in practice. 

In order to see what happens near the lower boundary of the conformal window, we 
need to define the gauge beta function in the non-perturbative region. Here we shall apply 
the ERG formalism. Because the RG equations can be obtained without recourse to the 
perturbative expansion, although some approximation is inevitable to solve them. The 
Wilson RG describes the flows in the infinite dimensional coupling space by integrating 
out the higher momentum modes scale by scale 6j. For a renormalizable theory, the flows 
converge towards a line called the renormalized trajectory (RT) in the continuum limit. The 
RT corresponds to the renormalized theory. The beta function for the renormalized theory is 
nothing but the scale transformation on the RT, and therefore is given in a non-perturbative 
way in principle [8[. 

Here it should be noted that the Wilson RG offers us a very suitable framework in seeking 
for the fixed points as well as the phase structure. It is also a great advantage to find out 
the anomalous dimensions at the fixed points directly from the RG equations j^]. It is also 
remarkable that the Wilson RG is applicable not only in the broken phase but also in the 
symmetric phase. Contrary to these features, the DS equations cannot be applied in the 
symmetric phase. Therefore the DS equations cannot examine the IR dynamics around the 
IR fixed point, while it is useful to evaluate the order parameters in the broken phases 



In the explicit calculation, we will incorporate the general four-fermi operators 12 



JJaL 



allowed by the symmetries into the effective Lagrangian, but truncate other higher dimen- 
sional operators. The essential point is to extend the RG equation for the gauge coupling so 
as to include the higher order corrections via the four-fermi operators. More explicitly, we 
give the RG flow equation for the gauge coupling by adding the gauge invariant corrections 
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via the effective four-fermi interactions to the two-loop perturbative beta function. This de- 
pendence on the four-fermi coupling enables us to take non-perturbative corrections in the 
gauge beta function. We also discard all the gauge non-invariant corrections as a primitive 
analysis. 

Thus we can define the gauge beta function containing non-perturbative corrections 
through non-renormalizable effective operators. The RG flow connecting the UV fixed 
point and the IR fixed point is found to be the RT representing the continuum limit of 
the asymptotically safe theory. Therefore the non-perturbative gauge beta function also 
possesses a UV fixed point. Meanwhile it will be also shown that this RT is not achieved by 
perturbative expansion of the RG flow equations. Furthermore The UV fixed point in the 
non-perturbative gauge beta function is found to merge with the IR, or the BZ, fixed point 
at the boundary of the conformal window. 

Recently several ansatz for the non-perturbative gauge beta function have been proposed 
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22] . In Ref. 22|, the fixed point merger is imposed as one of the ansatz. Contrary to 



these approaches, we do not make an ansatz, but try to consider the continuum limit and 
the beta function in a tractable approximation scheme. Therefore we do not intend to seek 
for the analytic forms of the non-perturbative beta functions. 

This paper is organized as follows. In section 2, we will discuss a toy example of a scalar 
theory to show how we can define the non-perturbative beta function by applying the ERG 
framework. Then, in section 3 we derive the RG equations for the SU(N C ) gauge theory 
with Nf massless flavors in the above mentioned approximation scheme. In section 4, we 
will show the numerical results obtained by solving the RG flow equations. First the fixed 
point structure and the anomalous dimensions of fermion mass operator at the fixed points 
are shown. The perturbative analysis for the gauge beta function are also examined and 
discussed. The non-perturbative gauge beta functions and the fixed point merger are shown 
in section 5. The last section is devoted to the conclusions and discussions. 



II. NON-PERTURBATIVE BETA FUNCTIONS GIVEN BY THE WILSON RG 

In this section, we would like to consider how to define the beta function in a non- 
perturbative way by using the ERG equations. We may define the Wilsonian effective action 
by integrating out the higher frequency modes by introducing momentum cutoff. Then all 
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operators allowed by the symmetry are induced in the effective action irrespectively of the 
renormalizablity. Then the effective action is represented by a point in an infinite dimen- 
sional parameter space (theory space). The Wilson RG follows response of the Wilsonian 
effective action under scale transformation and, therefore, defines a RG flow for each theory 
in the theory space. 

In general, the flows converge towards a special flow called the renormalized trajectory 
(RT) in the infinite cutoff limit (continuum limit). The renormalizable theories are described 
in terms of the renormalized trajectories. Accordingly the beta functions of the renormalized 
parameters are given by the scale transformation on the RT in the Wilson RG framework. 
It is noted that the Wilson RG equations may be given without perturabative expansion, 
although approximations such as operator truncation must be performed. Therefore the 
beta function defined on the RT contains non-perturbative informations, partly though. In 
the followings, let us consider these aspects explicitly by examining the Wilson RG equations 
for a massless scalar 4 theory. 

We truncate the Wilsonian effective action as 

**±iw> - ^ - (i) 

where we discarded the mass parameter just for simplicity. Then the theory space is two 
dimensional and is spanned by (A4, Xe). Here we may use the Wetterich equation j^] given 
in terms of t = In A by 

d t r A [0] = ^STrd t R A (rf [0] + RaY' , (2) 

where stands for the cutoff effective action and T A denotes the second functional 

derivative with respect to the scalar field <fi. The function R\{p) is the regulator function, 
which plays a role of cutoff around p 2 ~ A 2 . Then the RG flow equations are found to be 8| 

A^f = aX\ - bX 6 , (3) 
A^ = 2A 6 - cX\ + 2dA 4 A 6 , (4) 

where the coefficients are dependent on the choice of the cutoff functions R^. If we take the 
sharp cutoff limit, then the coefficients are given as 

a = 9A, b= 10A, c = 27 'A, d = yA, (a = JL) . (5) 



In Fig. 1, the RG flows are shown by the dotted lines. It is seen that the flows converge 
rapidly towards the red line, which shows the RT obtained in the continuum limit js]. Once 
the RT is found, then irrelevant coupling Xq is given as a function in terms of the renormalized 
coupling A4, or Xq = Ag(A 4 ). Therefore the beta function is given explicitly as 



At(A 4 ) = a\ 2 4 - 6Aq(A 4 ). 



(6) 
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FIG. 1: RG flows obtained by solving Eq. ([3]) and dH) are shown by the dotted lines. The dashed 
lines represent the perturbative RTs obtained in each order. The red line (RT.) is the non- 
perturbative RT obtained numerically in the continuum limit. 

Here we would stress that this beta function contains non-perturabative corrections 
through the non-renormalizable coupling X e . In order to see that, we compare this beta 
function with those obtained by perturbative renormalization. If we solve the RG equations 
by expanding the couplings Aj(z = 4, 6) as 

A, = Af + AA« + A'xf + • • ■ , (7) 

then the solutions are easily found order by order. 
At the level of O(A ), the solut ions are found to be 

Ai 0) (A) = Af(Ao), (8) 
Af(A) = (A/Ao) 2 Af(A ). (9) 
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Here A4 (A) is fixed by imposing the renormalization condition, while Ag°^(A) = in the 
continuum limit, A — > 00. 

Similarly we may proceed to 0(A 1 ). The RG equations 

A^ = 2A( 1) -cAi 0)3 + 2rfAi 0) A(° ) , (10) 



give the solutions 



Ai X) (A) = alnA A f (A ) 2 + A ( 1 )(A ) ; (n) 



4"(A) = \ 



1 



A r 



Ar(Ao) + ( Tq ) A^(A ). (12) 



The logarithmic divergence in Eq. ( ITT]) is removed by the 1-loop counter term A^^Aq). In 
the continuum limit, it is seen by noticing Eq. (Q that 

A?>(A) |Af(A) 3 . (13) 

This is nothing but the relation on the RT in the first level of perturbation. In order to seek 
for solutions in the higher orders of A , we may iterate the similar calculations. Then it is 
found that the RT is given as 

A* = < lXl + ^(3a-2d)Xi- < ^(-12a 2 + 3bc+Uad-4d 2 )Xl 
2 4 8 

+ -(30a 3 - 18abc - 47a 2 d + lObcd + 24ad 2 - 4d 3 )\l + ■■■ . (14) 
8 

In Fig. 1, the dashed lines represent the perturbative RTs obtained in each order. The 
red line (R.T.) is non-perturbative RT obtained numerically in the continuum limit. It 
seems that the perturbative RTs approach to the non-perturbative RT obtained by solving 
RG flow equations directly. Thus the RT given by solving the ERG equations is found to 
contain non-perturbative series. 

In Fig. 2, the non-perturbative beta function defined by Eq. is given by the red line 
(NP /3 4 ). The dashed lines represent the perturbative beta functions obtained in each order. 
It is seen that these beta functions converge to the non-perturbative beta function. Thus 
the beta function containing non-perturabative informations can be easily found by solving 
the Wilson RG equations. 




FIG. 2: The non-perturbative beta function defined by Eq. © is given by the red line (NP ^4). 
The dashed lines represent the perturbative beta functions obtained in each order. 

III. WILSON RG FOR THE SU(N C ) QCD WITH N f FLAVORS 

We consider the SU(N C ) gauge theory with Nf vector-like massless fermions. The theory 
satisfies the global chiral symmetry SU(Nf)L x SU(Nf)n as well as the U{l)v and the 
parity symmetry interchanging the left-handed and the right-handed fermions. Therefore 
the Wilsonian effective action I\ also should maintain these symmetries in addition to the 
SU(N C ) gauge symmetry. 

It has been found that the phase structure of chiral symmetry is revealed by incorporating 
the four-fermi operators in the effective Lagrangian. Traditionally the chiral symmetry 
breaking has been studied by means of the DS equations. However, The Wilson RG has the 
following advantageous features to the DS approach. First it is straightforward to find out 
the phase boundary as well as the fixed points. Second, it is rather easy to improve the level 
of approximation by including more operators in the Wilsonian effective action. In practice, 
it has been shown that the non-ladder corrections can be treated equally with the ladder 
corrections in QED and the gauge dependence, which bothers us in the DS approaches, 
is remarkably improved Third, it is naturally performed to incorporate running effect 
of the gauge coupling in the Wilson RG framework, while one may apply the Higashijima 
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approximation in the DS equations 23]. 



Therefore, it is important to include four-fermi operators in order to examine the chiral 
phase structure and the conformal window of the SU(N C ) gauge theory. In this paper we 
take only the four-fermi operators in the Wilsonian effective action Ta, but discard all other 
higher dimensional operators as the first study; 

r A = J d 4 x (--Lf*f a v + i> f ilp^ + £ 4/ ) , (15) 

where the color index a of massless fermions represented as ip a f,(a = 1, ■ ■ • ,N c ,f = 
1, • • • , Nf) is suppressed and £ 4 / denotes the four-fermi operators. Strictly speaking, gauge 
variant operators are also induced through the RG evolution due to momentum cutoff. 
Therefore we should include these operators in the effective action as well. However, we 
choose an easy way by discarding all gauge variant corrections. We shall come back this 
point later. 

The four-fermi operators must be invariant under the SU(N C ) gauge transformation, 
the SU(Nf)i x SU(Nf)n chiral transformation and the parity. We may write down vari- 
ous invariant operators, however some of them are transformed to each other by the Fierz 
;ransformation. It is shown in Appendix that there are 4 independent four-fermi operators 
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20| . which are given, for example, by 35] 



where 

o s = -iL.wnj: = i - ^75^'^75^] , (17) 

Oy = I.r"I. j I, r! ,I' + {L^R)= l - [jvfWjl^ + ^W^iT^*] , (18) 

Ovi = 2/.r / 7.7,y / ,/," = ~ [(^yy) 2 - $i7V) 2 ] > (19) 

V2 = ila^Uf + {L^R) = l - p T V) 2 + (^7 M 75#) 2 ] • (20) 

Here the color indices of fermions are always contracted and i,j denotes the flavor indices. 
The left(right)-handed fermions are represented by L m = ip™ (R m = if) c } 

The first operator Os plays an essential role for chiral symmetry breaking [5]. The 
coupling Gs grows rapidly in the broken phase and diverges, which indicates chiral symmetry 
breaking of 

= M 3 5{. (21) 



This non-trivial order parameter is expected just like in the mean field approximation of the 
Nambu-Jona Lasinio (NJL) model, and then the flavor symmetry SU(Nf)i, x SU(Nf)n is 
broken to SU(N f ) v . 

Now we shall derive the ERG equations for this system. By expanding the Wetterich 
equation by operators, we may obtain the RG equations for various effective couplings. The 
trace form of the quantum corrections in the ERG equation given by Eq. implies that 
they are just one-loop ones, although the vertices include non-renormalizable interactions. 
Therefore we may calculate the RG equations for the effective couplings by evaluating various 
one-loop diagrams with shell integration for the loop momentum. 

For the four-fermi couplings, the one-loop diagrams are represented as in Fig. 3, where 
four-fermi vertices represent any of the operators given by Eqs. ( lTTll2"Uj) . The wavy line 
stands for the gauge propagator. Thus we need to calculate rather many diagrams. So 
we shall summarize details of the calculations in Appendix and present only the results here. 



FIG. 3: The RG corrections for the effective 4-fermi couplings are illustrated diagrammatically. The 
bold and wavy lines represent propagators of the chiral fermions and the gauge bosons respectively. 

We use the Landau gauge propagator and evaluate the loop integrations in the sharp 
cutoff limit. We also use g { = Gi/ATT 2 (i = S,V,V1,V2) and a g = g 2 /(A7c) 2 . Then the 
resultant RG equations are found to be 




A 



dgs 



2g s - 2N c g 2 s + 2N f g s g v + §9s9vi + ^9s9v2 



dA 




(22) 



A 



dg v 



2g v + (N f /4)g 2 + {N c + N f )g 2 v - Qg v g V2 



dA 




(23) 



A 



dgyi 
dA 



2sVi - (1/4)0! - 9s9v - 39vi - N f9s9v2 + 2(N C + N f )g v g vl 




(24) 
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A ^T = 29V2 ~ ^ ~ NcN ^ + ^ NcN f ~ 2)9 ^ ~ N f9s9vi 

+2{N c N f + l)g v g V2 + 6(g v + g V2 )a g - ^ (3 + -^J a% (25) 

where 02(F) denotes the quadratic Casimir of color reprerentation of the fermions and 
is given explicitly by (N 2 — 1)/2N C . In deriving these equations we do not make further 
approximations, and therefore valid for any N c and Nf. 

We need to solve these differential equations coupled with the RG equation for the gauge 
coupling. Therefore the RG flows are given in the five dimensional parameter space and 
the flow diagram becomes rather complicated. Since our present purpose is qualitative 
understanding of the non-perturbative RT and the beta function, let us coincide the RG 
equations in the large N c and Nf limit. This limit is taken by rescaling as 

N c gs(V) -»■ 9s(V), (26) 
^c9vi(V2) ->■ gvi(V2), (27) 
N c a g a g , (28) 

with keeping the ratio r = Nf/N c . Then it is seen that the first two equations, (1221) and 
( 123]) are reduced to 
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A 1K = 2gs ~ 2g2s + 2rgsgv ~ Qgsag ~ (29) 

= *9v + - A 9l + (1 + r)g v - 3 -a 2 . (30) 



It is noted that four-fermi couplings gy\ and gv2 decouple from the above equations. There- 
fore we may solve only three equations for the couplings gs,9v and a g in the large N c and 
N f limit. 

Next we consider the RG equation for the gauge coupling. The 2-loop beta function for 
the SU(N C ) QCD is given by 

da g 
dp, 

where a g = g 2 /(47i) 2 and 



= ^—1 = -2b a 2 g - 2ha% (31) 



bo = yiV c - |jV>, (32) 
34 „ fN 2 -l 10 



k = -jN; - N, + T N C ). (33) 
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FIG. 4: The lowest order correction to the gauge vertex depending on the effective four-fermi 
interactions. 



Here we intend to incorporate non-perturbative corrections into the gauge beta function 
in the same way performed for the scalar field theory in section 2. It is thought that it 
is important to find dependence on the four-fermi couplings in the RG flow equation for 
the gauge coupling. Then some non-perturbative corrections are taken in the gauge beta 
function through the effective four-fermi couplings. 

However, it is not straightforward to apply the RG flow equations for the gauge theories. 



One problem is that the Wilson RG does not respect manifest gauge invariance 36[ . Indeed, 



we may make use of the modified Slavnov- Taylor identities in order to derive the RG flow 
equations for gauge invariant theories [25]. Another problem is that we need to deal with 
fairly complicated flow equations involving effective operators with higher derivatives in 
order to reproduce even the two-loop beta function. So we give up to apply the Wetterich 
equation for the gauge coupling directly and, in stead, derive the RG flow equation along 
the following considerations. 

First, we substitute the part of the flow equation depending only on the gauge coupling 
by the two-loop perturbative beta function given in Eq. f l31~|) . Then we add corrections 
depending on the four-fermi couplings as the non-perturabative part of the beta function. 

We may think of one-loop corrections shown in Fig. H] inserting the four-fermi vertices 
as the non-perturbative corrections. The explicit calculation of the shell mode integration 
leads to the extra contribution to the gauge beta function, 

5a g = -r(g s - 2g v )a g 5lnA, (34) 

in the large N c and Nf limit. However, it is noted that these corrections should be forbidden 
in the Landau gauge calculation, if the gauge invariance is protected. It may be also con- 
firmed by the Slavnov- Taylor identity that the vertex correction is indeed gauge variant jl9| . 
Therefore, we discard these corrections from our RG flow equation for the gauge coupling. 
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FIG. 5: The two-loop correction to the vacuum polarization of the gauge field depending on the 
effective four-fermi interactions. 



However, we expect that there should exist gauge invariant corrections to the gauge beta 
function through the effective four-fermi couplings. So let us consider the QED case. The 
higher order corrections to the gauge beta function come from the vacuum polarization 

ri 

diagrams. Especially we may focus on the ladder type diagrams shown in Fig. [5] [37]. In the 
Wilson RG picture, a part of such contributions can be taken in by replacing the middle 
box diagram with the effective four-fermi operators as shown in Fig. |5j 

The box diagrams with fermions of the same chirality induce the effective four-fermi 
operator given by 



efl 



A 2 
2GV 

A 2 



[L n »T A U L jltl T A U + (L^ R)] 



(35) 



where A denotes the scale of the effective theory and Gy is the induced four-fermi cou- 
pling. In the second line we performed the Fierz transformation (see Appendix). Then 
the four-fermi operator may be regarded as the current- current interaction connected by a 
gauge propagator with momentum scale A. Therefore such higher order contributions in the 
vacuum polarization diagrams may be sum up by replacing the gauge boson exchange in 
the two-loop diagram by (g 2 + 2Gy)D tlv , where D^ v denotes the gauge boson propagator. 
The two-loop vacuum polarization is calculated with dimensional regularization in the gauge 
invariant way. Thus the vacuum polarization li^ including the effective interaction is found 
out to be 



IW<z) = —Nf 



4 



-a g + 2a + a g g v 



3 " " 

where the epsilon pole in the dimensional regularization is replaced by In A. 



(36) 
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Consequently we may derive the RG flow equation for the gauge coupling a g as 

A ^A = ^ + 2rgW v (37) 
where the rescaling for the large N c and Nf limit has been performed and the two-loop beta 
function (3 g is given by 

/f = -|(H-2r)aJ-|(34-13r)aJ, (38) 

in this limit. It is noted that the second term of Eq. (1371) corresponds to — bX^ in Eq. (J6]) 
defined for the scalar toy model, which contains non-perturbative informations. We used 
the two- loop beta function p g as the perturbative part, since the non-perturbative part 
contains corrections only of more than the two-loop order. We may apply the three- or 
four-loop beta function 0, Q to the perturbative part as well, although the corrections 
may be overlapping with the second term partly 38(. In this paper we perform explicit 
calculations of the RG flow equation given by Eq. (1371) only. 



IV. RG FLOWS AND THE RENORMALIZED TRAJECTORIES 

In this section we show the RG flows obtained by solving numerically the RG flow equa- 
tions in the large N c and Nf limit given by f )29|) . f l30]) and fl37|) . Before studying the gauge 
theories, let us examine the flows of the pure four-fermi couplings by setting a g = 0. The 
RG flows in the gs~gv plane are shown in Fig. |6j It is seen that the RG flows are separated 
into two parts by the phase boundary. The chiral symmetry is broken in the upper region, 
since gs diverges towards the infrared direction. On the phase boundary there is a UV fixed 
point, where gy takes a negative value. The IR fixed point is trivial. 

At the UV fixed point the four-fermi operators are recombined into a relevant one and 
an irrelevant one. The relevant operator, 

O rel = V Os + O v , (39) 

where the coefficient is i] = 2 (—2 — \Jr 2 — r + 4)/r. The scaling dimension of this operator is 
found to be just 2 at the fixed point. Though the relevant direction is given by a straight line 
as seen in Fig. [6j the line will be found to be curved in the presence of the gauge interaction. 

Next we examine the RG flows of the gauge theories. It is immediate to find the fixed 
points, since they are just given by zero points of the RG flow equations f l29|) . fl30|) and f )37|) 
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FIG. 6: Aspect of the RG flows of the pure four-fermi interactions in the gs-gv plane. 



39[. In Fig. [7J the fixed points are shown in the 3 dimensional coupling space for various 
numbers of r = Nf/N c . In the conformal window the IR fixed point carries non-trivial 
couplings and a UV fixed point also appears. As reducing the flavor number, these fixed 
point couplings become strong. Eventually the UV fixed point and the IR fixed point merge 
with each other and disappear around r ~ 4.05. 
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0.10 
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FIG. 7: UV and IR fixed points for r = N f /N c = 16/3, 15/3, 14/3, 13/3 and 4.05 



The scaling dimensions of the fermion operator at the vicinity of the IR and the UV fixed 
points are important physical quantities to be measured. Especially, the eigenvalues of the 
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RG flow equations linearized around the fixed point are nothing but the scaling dimensions 
of the four-fermi operators. We are interested in the most relevant dimension, since it shows 
the (ir) relevance of the four-fermi operators at the fixed point. Fig. [H] shows the dimensions 
evaluated at the UV (blue) and the IR (red) fixed points in the cases of various flavor 



numbers 
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131 ] . The scaling dimension obtained in the large N and ladder approximation 



is also shown for comparison. It is seen that the dimensions are almost same as those 
evaluated in the ladder approximation. It is noted that the dimension should be 4 exactly 
at the edge of the conformal window, since the four-fermi operator becomes marginal by 
fixed point merger there. 

In the present scheme, the anomalous dimension of fermion mass is simply represented by 
the diagrams given in Fig. js), [jjj]. The crossed vertex stands for the mass insertion. It is 
immediate to evaluate these corrections at the fixed points and the scaling dimension of the 
fermion mass operator A^, are shown in Fig. [10] for various flavor numbers in the conformal 
window. We also present the results obtained in the large N and ladder approximation for 
comparison. It is found that the scaling dimensions at the UV fixed points differ from the 
ladder results significantly, while there is not much differences at the IR fixed points jl2]. It 
is noted that sum of the anomalous dimensions at the IR fixed point and the UV fixed point 
is more than the value, —2, which is expected in the large N c analysis 28] . This is because 
the corrections taken in the RG flow equations contains not only planer diagrams but also 
non-planer ones due to large Nf contributions, which influence the anomalous dimension 
significantly. 

Next we consider change in behavior of the RG flows near the lower edge of the conformal 
window. In Fig. [11] and Fig. [12] the RG flows of SU (3) gauge theories are shown in the 
case with Nf = 13 and Nf = 12 respectively. The fixed points are shown by red points. For 
Nf = 13, the theory lies within the conformal window and the UV and the IR fixed points 
with non-trivial gauge couplings exist. Meanwhile the theory with Nf = 12 is slightly out 
of the conformal window and these fixed points have annihilated in pair. 

First let us discuss the RG flows of the theories in the conformal window. As seen 
in Fig. [11] there is the phase boundary and the chiral symmetry is broken in the upper 
region. The flows in the lower symmetric region run into the IR fixed point. Then the flows 
converge towards a line connecting three fixed points; the trivial fixed point, the IR fixed 
point and the UV fixed point. This line is nothing but the renormalized trajectory. The 
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large N c and ladder 



i - 

4.2 4.4 4.6 4.8 5.0 5.2 5.4 

r = N;/N c 

FIG. 8: Scaling dimensions of the most relevant four-fermion operator at the IR and the UV fixed 
points in the conformal window. 




FIG. 9: Diagrammatic representation for the anomalous dimension of fermion mass operator. The 
crossed vertex stands for the mass insertion. 



N f = 13 



N f = 14 1 



N f = 15 



N f = 16 




large N c and ladder 



4.2 4.4 4.6 4.8 5.0 5.2 5.4 



r = Nj/N c 



FIG. 10: Scaling dimensions of the fermion mass operator at the IR and the UV fixed points in 
the conformal window. 
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9v 




FIG. 11: Aspect of the RG flows in the 3 dimensional theory space in the case of Nf = 13. The 
fixed points are marked by red points. 



9v 




FIG. 12: Aspect of the RG flows in the case of Nf = 12. 

flow starting from the UV fixed point for the IR one corresponds to an asymptotically safe 
theory, while the flow approaching the IR fixed point from the trivial fixed point corresponds 
to an asymptotically free theory. The both theories can be defined as the continuum limit 
and, therefore, renormalizable. 

On the other hand, it is seen in Fig. [12] that the non-trivial fixed points and the phase 
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boundary disappear just below the conformal window. The four-fermi coupling g$ diverges 
in all flows, therefore the theories in the entire region belong to the broken phase. The 



RG flows converge into a RT corresponding to the the asymptotically free theory 40J. The 
RT is unique and extended to infinity. When the flavor number Nf is lowered further, the 
RG flows show the basically same structure and the RT remains in the QCD with a small 
number of flavors. Now it is clearly seen from these two figures that the RT connecting the 
trivial fixed point and the UV fixed point passing through the IR fixed point for the theories 
in the conformal window is transformed continuously to the asymptotically free RT below 
the conformal window. 

However, when we look into the RG flows in the conformal window more carefully, then 
these show a somewhat curious feature as follows. The RG flows obtained in the case of 
Nf = 15, which is near upper bound of the conformal window of the SU(3) gauge theory, 
is shown in Fig. [13j The flows are projected on the (a g ,gs) plane. There is the UV fixed 
point as well as the IR fixed point and, therefore, the flow connecting them is a RT, while 
the flow connecting the trivial fixed point and the IR fixed point is also another RT. 

We note first that the RG flows in the symmetric phase converge into a distinct line from 
the RT, when the gauge coupling approach the IR fixed point from the strong side. The 
RT of the asymptotically free theory is unique. Therefore, there seem to exist two different 
theories remaining in the continuum limit. 

When the IR fixed point couplings are sufficiently weak, then the perturbative analysis 
is reliable around the fixed point. So we examine the RT by solving the RG flow equations 
given in (12 9j) . ( 1301 and ( 1371) by perturbative expansion. We may carry out the similar analysis 
performed for the scalar model and find the following relations in the continuum limit, 



9s = -^l - ^(-3 + 2&oK + ^(90 - 120&o + 486jj - - 3r)a g + ■ • • , (40) 
9v = \°? g - h a 3 g + JL(-3 + 96b 2 - 32&i - 30r)aJ + ■■■. (41) 



In Fig. [T3], the RTs obtained in this manner up to the 8-th order are also shown order by 
order. It is seen that the asymptotically free RT is perfectly described by the perturbative 
solutions. However the series of these RTs converge to a line distinct from the RT connecting 
the UV fixed point beyond the IR fixed point. This shows that the theory described by the 
RT connecting the UV fixed point cannot be obtained by perturbative renormalization. 
Therefore we may say that the asymptotically safe theories are truly non-perturbative in 
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this sense. 

While the perturbative RTs seem to converge around the IR fixed point, they do not seem 
to converge for a larger gauge couplings. Therefore it may be expected that the perturbative 
RT does not represent the true continuum limit in the non-perturbative sense. When the 
flavor number exceeds the conformal window, Nf > (ll/2)iV c , the theory becomes IR free 
and behaves like QED The continuum limit of 4 dimensional QED has been discussed 



in the context of ERG and absence of the continuum limit is pointed out 
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FIG. 13: The RG flows and the perturbative RTs obtained in various order calculations are shown 
in the a g -gs plane for Nf = 15. The asymptotically safe RT is shown by a bold line. 

As the flavor number Nf is reduced and the fixed point couplings become strong, the 
perturbative RT also becomes close to the non-perturbative one. In Fig. and in Fig. [151 
the projected RG flows are shown in the case of Nf = 13 and Nf = 12 the for SU(3) QCD 
respectively. There the RTs obtained by the perturbative expansion are also shown by the 
dashed lines order by order. Still the asymptotically free RT is found to be described well 



by the perturbative series 42]. 



It is seen in the figures that the perturbative RT gets closer to the RT passing through 
the UV fixed point for Nf = 13 and becomes indistinguishable from the asymptotically 
free RT for Nf = 12. So it is supposed that the perturbative RT and the non-perturbative 
RT become identical with each other out of the conformal window. These statements are 
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somewhat speculative, although we may naively expect them by looking at the figures. Of 
course, it is desirable to prove them mathematically, but it is beyond our present purpose. 



0.5 r 




FIG. 14: RG flows and perturbative flows in the continuum limit projected on the a g -gs plane for 
AT, = 13 



V. NON-PERTURBATIVE BETA FUNCTIONS 

Now we discuss the non-perturbative gauge beta functions for the many flavor QCD in 
and out of the conformal window. As discussed in section 2, the non-perturbative beta 
function can be given by the scale transformation on the RT. It was also shown that the 
dependence on the non-renormalizable couplings in the sense of the perturbative expansion 
is essential to obtain the non-perturbative beta function. Now the RG flow equation for the 
gauge coupling (1371) contains a term dependent on the four-fermi coupling g v , which reflects 
the non-perturbative corrections. 

First we examine the beta functions for the SU(3) gauge theory with Nf = 15, which 
lies in the conformal window. As we discussed, we can define the non-perturbative RT 
in addition to the perturbative RT, which seems to be distinct in the conformal window. 
The beta function depends on which RT we choose to define it. We may define the non- 
perturbative beta function given by Eq. (j57|) based on the non-perturbative RT, which 
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FIG. 15: RG flows and perturbative flows in the continuum limit projected on the a g -gs plane for 
N f = 12 

transforms to the RTs out of the conformal window smoothly. In Fig. [16] the beta functions 
are shown by the bold lines. The figure also presents how the non-perturbative beta function 
changes it's form as the flavor number is reduced. It is seen that the gauge beta function 
possesses a UV fixed point in addition to the IR fixed point in the conformal window and 
that these fixed points merge with each other and disappear at the edge of the conformal 
window. Such a behavior of the beta function has been discussed in more general contexts 
by Kaplan et. al. recently [n| • Our analysis shows that the "conformality lost" may realize 
in the gauge beta function of QCD with many flavors. 

On the other hand, we may also define the perturbative beta functions by substituting the 
solution given by Eq. ( l4Tj) into the non-perturbative beta function defined by Eq. (1371) . The 
beta functions are shown by the dashed lines in Fig. [T6j It is seen that the perturbative beta 
functions are almost same as the two-loop beta functions, which are shown by the dotted 
lines. Especially it is noted that no sign of the UV fixed point is seen for the perturbative 
beta functions. Although the RT obtained by the perturbative calculation approach the 
non-perturbative RT near the edge of the conformal window, the beta functions are clearly 
different. Thus we find that appearance of the UV fixed point in the gauge beta function is 
purely due to non-perturbative effect. 
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FIG. 16: The non-perturbative gauge beta beta functions are shown by the bold lines in some 
cases of different flavor numbers. The red points represent the fixed points. The perturbative beta 
functions and the two-loop beta functions are also shown by dashed lines and the dotted lines 
respectively for comparison. 



One may wonder that disappearance of the fixed point in the gauge beta function, or 
lost of the scale invariance, is due to fermion mass generation in the chiral symmetry broken 
phase. Furthermore, it may be thought that the bending behavior of the gauge beta function 
out of the conformal window is induced by fermion decoupling due to mass generation. How- 
ever, we may also make the following objections. Suppose that the beta function possesses 
only the IR fixed point and the gauge coupling a gm is slightly less than the critical gauge 
coupling a 9cr for the chiral symmetry breaking. According to the above picture, the beta 
function should bend in the region of a g > a 9ci by fermion decoupling. Then a UV fixed 
point appears in the beta function inevitably. However this is forbidden in the broken phase. 
The gauge couplings at the fixed points a 9IR uv should be less than the critical coupling a 9cr 
and this is exactly the case realized in our analysis. Obviously the bending behavior of the 
non-perturbative beta function between the IR fixed point and the UV fixed point cannot 
be induced by fermion decoupling, since the theory stays in the symmetric phase. 

In our RG anlaysis chiral symmetry breaking occurs when the four-fermi couplings goes 
to infinity. Therefore the scale of dynamically generated fermion mass should be much less 
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than the scale range of the RG flows in our analysis. Accordingly the bending behavior 
of the non-perturbative beta function out of the conformal window is not caused by chiral 
symmetry breaking either. It is an interesting problem to evaluate the dynamical scale in 
our RG framework, though it is beyond our present scope. Then we may explicitly discuss 
the gauge beta function with fermion decoupling effect and clarify the differences from the 
non-perturbative beta function obtained in this paper. 

VI. CONCLUSIONS AND DISCUSSIONS 

In this paper we considered the gauge beta functions of the SU (N c ) gauge theories in the 
conformal window. So far, dynamics of the chiral symmetry breaking near the conformal 
window has been investigated by means not only of the DS equations but also of the lattice 
MC simulations. The phase structure and the IR/UV fixed points in the conformal window 
have been already examined by the ERG approach as well. So far the analytic approaches 
adopted the perturbative gauge beta functions. This paper presents the first study in order to 
understand the non-perturbative aspect of the gauge beta function in and near the conformal 
window. 

It is shown explicitly that the non-perturbative beta function can be evaluated without 
perturbative expansion in the Wilson RG. Here, of course, non-perturbative does not mean 
exact, since we cannot help but performing some approximations in the practical analysis 
of the RG flow equations. We incorporated 4 independent four-fermi operators allowed by 
the symmetry in the Wilsonian effective action, but discarded all other higher- dimensional 
operators. Then the RG flow equations for these four-fermi couplings are derived explicitly 
in the sharp cut-off scheme. In practice, we analyzed the flow equations in the large N c and 
Nf limit, in order to discuss qualitative behavior of the gauge beta function. 

We extended the RG flow equation for the gauge coupling so as to include the higher 
order corrections through the four-fermi operators. Explicitly we used two-loop beta function 
as the perturbative part and added the non-perturbative corrections through the effective 
four-fermi coupling. This is the most important point in our analyses. We discarded gauge 
variant corrections, which are induced due to momentum cutoff, just to make the analysis 
simple. Among all we neglected the vertex corrections with the four-fermi operators. 

We obtained the RG flows in the 3 dimensional coupling space by solving the flow equa- 
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tions numerically. In the conformal window, there exists a UV fixed point on the critical 
surface of the chiral symmetry in addition to the IR fixed point. At the edge of the confor- 
mal window, these fixed points merge with each other and disappear. Apparently the flow 
connecting the UV fixed point and the IR fixed point is a RT, which describes an asymp- 
totically safe theory. This RT is found to be transformed continuously to the RTs of the 
asymptotically free theories out of the conformal window. 

On the other hand we may take the continuum limit of the perturbative solutions of the 
RG flow equations. However it is found that the perturbative RTs do not converge to the 
asymptotically safe RT in the conformal window, although convergence to the asymptotically 
free theory is remarkably good. In other words, perturbative renormalization cannot describe 
the asymptotically safe theory. Moreover the series of the perturbative RTs itself does not 
seem to converge at the strongly coupled region. It is seen also that the perturbative RT 
coincides with the RT obtained by RG flows in the theories below the conformal window. 

We defined the non-perturbative gauge beta function, which gives the scale transforma- 
tion on the RT. Then it is found that the gauge beta function possesses not only the IR 
(BZ) fixed point but also a UV fixed point and these fixed points merge with each other at 
the edge of the conformal window. Meanwhile the beta function evaluated by the pertur- 
bative analysis of the RG flow equations does not show the UV fixed point. In addition, 
the perturbative beta function keeps an IR fixed point even below the conformal window. 
Thus it is said that the phenomenon of fixed point merger is owing to the non-perturbative 
corrections taken in the gauge beta function. 

Recently Kaplan et. al argued that the BKT type phase transition occurs in such situa- 
tions in general . For gauge theories in the conformal window, this phase transition has 



been known as the Miransky scaling 



18l | . In this paper we gave the explicit beta functions 



showing the "conformality lost" in the SU(N) gauge theories. 

The order parameters of chiral symmetry breaking can be evaluated by using the RG 
flow equations in the broken phase. It would be an interesting problem to verify the scaling 
properties with respect to the flavor numbers and compare with the BKT type scaling 
expected by the gauge beta function j^J. Further studies along this line of consideration 
will be reported elsewhere. 

Lastly let us make a comment on the beta function for the N — 1 supersymmetric QCD. 
It has been known that there is also the conformal window in the supersymmetric QCD and 
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the region is exactly given by (3/2) N c < Nf < 3N C [32] . It is remarkable that the anomalous 
dimension of the mass parameter at the fixed point is related with the charge of R-symmetry 
carried by the quark superfields and, therefore, is determined exactly. Accordingly the 
anomalous dimension at the fixed point is unique. This seems to imply that the fixed point 
is also unique, or that the UV fixed point is absent. 

Furthermore the N=l supersymmetric QCD is known to satisfy the so-called exact beta 
function |33| , which is given in terms of the anomalous dimension of the quark superfields. 
If the anomalous dimension is a monotonic function of the gauge coupling, then the beta 
function also has a unique fixed point. We cannot know how the gauge coupling at the fixed 
point changes in the conformal window, since we do not know the anomalous dimesion as a 
function of gauge coupling except for the perturbative calculation. However uniqueness of 
the fixed point suggests that the fixed point value moves to infinity in order to transform 
to the beta function without a fixed point continuously. If magnetic coupling of the dual 
gauge theory behaves like inverse of the gauge coupling, then the infinitely strong fixed 
point is thought to be allowed 17j. This feature seems quite contrasting with the non- 
perturbative beta function discussed in this paper . Thus extension of the present study 
to the supersymmetric cases will be also quite interesting and remains for future studies. 
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Appendix A: Invariant four-fermi operators 



In this section we show that the four-fermi operators given by ( 1T71I201) offer the inde- 
pendent basis of the invariant four-fermi operators 12 L The set of independent four-fermi 
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20[, however, the operators given 



operators has been also discussed previously in Ref. 
there are somewhat different from ours. So, we make a brief explanation on the invariant 
four-fermi operators and also present some useful formulae. 
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The effective four-fermi interactions of massless SU(N C ) gauge theories should satisfy the 
global SU(N/)l x SU(Nf)n chiral symmetry and the parity interchanging the left- and the 
right-handed fermions in addition to the SU(N C ) gauge symmetry. Here we use a, 5, • • • for 
the color indices and i,j, ■ ■ ■ for the flavor indices and represent components of a left- and a 
right-handed chiral fermion by = L a% and ip^ = R a% respectively Then it is found that 
any invariant operator is given by a combination of the following operators; 

2L ail »L ai R bjl ^R b i = X - [^ a a^ ai f ~ (Vrf75<n 2 ] , (Al) 
2L ai YL bi R bjl ,R a > = ~ ^ ai7 V%7^' - ^l^bilnlsV 3 ] , (A2) 
L ai ^L ai L bjlft L bj + (L^R) = l - [{^r'f + (^7,r 1 ) 2 ] , (A3) 
L ail »L u L bjllM L aj + (L <-> R) = X - [Aa^^jl^ + ^Yl^^l^r 3 ] , (A4) 

L ai YL a ^L bjlfl L bi + {l^r) = 1 - [^r^bji^™ + ^i^ aj ^i^ u ] , (A5) 
L ail ^L bjlfl L ai + (L R) = X - [^a^A^r* + ^i^^i^r 1 ] ■ (A6) 

Then we may obtain following formulae of the Fierz transformation, which are useful in 
calculations of the RG flow equations for the four-fermi couplings; 

^17^2^37^4 + V'l7 /i 75V'2^37A t 75^4 = ^ll^A^l '^2 + ^17^75^4^30 'pjsfa, ( A 7) 

and 

^I^A ~ ^175^2^375^4 = [^17^4^37/^2 ~ ^17^75^4^37^75^2] • (A8) 

By using these identities the operators given in (lAlHA6j) are rewritten into 



(Al) = - {^ b ^ b3 r i - ^75^'^75<T) = -a ai R b m b] L a \ (A9) 
(A2) = - ($art aj M hi -ipaijsr^bjjs^) = —4:L ai R aj R bj L bi , (A10) 
(A3) = ~ (i> aa ^ bj $ bj l^ ai ~ ^l^^l^r 1 ) = (A6), (All) 

(A4) = i (^r^bji^ - ^r^r^bji^ 1 ) = (as), (A12) 

Therefore it is seen that there are 4 independent operators. We have selected the operators 
in which the color indices are contracted with in the bi-linears among all. Then the operators 
given by (ll7H2(jp are chosen as the independent basis of the invariant 4-fermi operators. It 
is also noted that the operator Og, which is the NJL type four-fermi operator, is found to 
be suitable to explore the dynamics of chiral symmetry breaking. 
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The Fierz identities containing the color generator matrix T A are also given by using the 
completeness identity, 

dimG 1 

2 E ( T X ( TA )l = « - w m - (A13) 

A=l c 

Then the formulae are given as follows; 

2 L t T A ^V RjT^Ri = -O s - ^O vl (A14) 

A c 

J2 UT A YL l L^^I? + (L^R) = ho v - ^0 V2 (A15) 

.4 c 

UT A -fU LjT^U + (L^R) = \o V2 - ^-O v (A16) 



These identities are also useful in deriving the RG flow equations given by Eq. 



Appendix B: RG flow equations for the four-fermi couplings 

In this section, we show a sketch of derivation of the RG flow equations given by 
in the sharp cut off limit. Corrections in the Wetterich equation ()2]) are represented as one- 
loop diagrams with effective couplings. Therefore the following diagrammatical derivation 
would be the easiest. The momentum integration is restricted into an infmitesimally small 
shell region by the derivative of the cutoff function d t R\- In the sharp cut off limit we may 
use 

dp = J^L8Qp\-A)6t, (Bl) 

where t = In A, and p stands for the Euclidean momentum. 

We need to evaluate corrections for the effective four-fermi couplings, and, therefore, con- 
sider the one-loop diagrams with four external fermions. However there are many diagrams 
to be concerned, and we shall show only typical two kinds of diagrams among them here. 

First we consider the diagrams using two effective four-fermi operators Os- There are 
some diagrams depending on the ways of contraction. Here we consider only diagrams shown 
in Fig. 17. The correction represented by the diagram (a) is explicitly given by 

(a) = _(^y w „/vte^'i£^*L' 



A 2 J C J * \i> 2 
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L{ Rj Lb 




Li Rj Lk 
(a) (b) 

FIG. 17: The typical examples of diagrams for the RG corrections to the effective four-fermi 
couplings involving 2 Gs interactions. 

Therefore this correction induces variation of the effective coupling Gs by 

SG S = ~^G S St (B3) 
The correction given by the figure (b) is given by 

GO = U^fX f dp L^-L k L k %* + (L^R) 



2 V A 2 

N- 



16tt 2 s A 2 

Therefore this generates the correction to Gy given by 



L-GlSt -Oy. (B4) 



Nt 

SG V = —^G 2 S St (B5) 



We may repeat the similar calculations, though there are 29 one-loop diagrams involving 
two four-fermi couplings in total. 

Another typical four-fermi diagram is given by the gauge interactions as shown in Fig. 18, 
which induces a correction of order of g 4 . The gauge propagator in the Landau gauge is 

= J 2 (a- " *f) ■ <B6) 

The correction given by the diagram (c) is evaluated as 

f' _ ifr _ iS l . 

(<0 = <f E J * hT A ^ j-T B YL k R l T B r-^T A YR'D^D^ 
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Rj Rj 
(c) 

FIG. 18: The typical example of diagrams for the RG corrections to the effective four-fermi cou- 
plings involving 4 gauge interactions. 

where we used the identity given by f)A13|) . Therefore this correction induces the following 
variations, 

SGs = -Te( N '--B& St > (B8 > 
iG « = T2-k& St - (B9 > 

We may find 4 different diagrams of order of g A including the crossed ladder diagrams and 
evaluate them in a similar way. 

Moreover there are 12 diagrams proportional Gig 2 (i = S, V, VI, VI). Thus we may derive 
the RG flow equations for the effective four-fermi couplings given by Eq. ( |22~H25l) . It is noted 
that the 3-points or the 4-points interactions of the gauge bosons do not contribute to the 
1PI diagram with 4 external fermions in the one- loop order. 
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